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Abstract 

The transfer matrix method is often used to calculate the input impedance of woodwind instruments but it ignores the 
possible influence of the radiated sound from toneholes on other open holes. In this paper, a method is proposed to account 
for external tonehole interactions. It is found that the external tonehole interactions increase the amount of radiated energy, 
* ^ . reduce slightly the lower resonance frequencies, and modify significantly the response near and above the tonehole lattice 

cutoff frequency. 

<^ , The results of simulations with the Finite Element Method, as well as experimental measurements, are presented and 

compared to the calculation with the method presented in this paper, confirming that the external tonehole interactions play a 
significant role in woodwind instrument. 

PACS: 43.75 Ef, 4320 Rz 

Keywords: radiation, musical instruments, reed instruments, clarinet 



1 Introduction 



A method to accurately and efficiently estimate the input impedance and resonance frequencies of woodwind instruments is 
of primary importance. The Transfer-Matrix Method (TMM) is typically used for this purpose (see e.g. Plitnik and Strong 



OlSll . Causse et al. yfl, Keefe ||l2|] ). because of its simplicity and efficiency. This method ignores internal interactions due to 
the coupling between the evanescent modes of nearby discontinuities as well as external interactions, which exists because 
the radiation impedance of each open tonehole is influenced by the radiation of sound from other toneholes. The problem 
of the response of woodwind instruments with external tonehole interactions was stated in a complete form by Leppington 



[15] and a method of solution was proposed by Kergomard lll3[l . using the mutual radiation impedance proposed by Pritchard 



191]. This method is based on the TMM for internal propagation with modified open tonehole radiation impedances to account 
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for interactions (referred to as TMMI). Preliminary experimental results were obtained by Springer 12011 for the case of holes 
spaced far apart in a pipe. No other validation of the method has been proposed. 

In this paper, we investigate the effect of external tonehole interactions in woodwind instruments with the TMMI. The first 
goal is to determine the validity of the method by comparing the results of TMMI calculations with Finite Element Method 
(FEM) simulations and wth measurements. We also compare these results with TMM calculations to show that most of the 
discrepancies are explained by the external interactions. 

The second goal is to apply the proposed TMMI method to the case of woodwind instruments in order to determine the 
importance of the effect on their acoustical properties and to judge whether or not it is necessary to account for those interactions 
when calculating the input impedance of woodwind instruments for design purposes. 

The theory of the TMM and TMMI, as well as the details of the FEM, are reviewed in the next section. This is followed by 
the presentation of the results for a tube with a regular array of holes (Sec. O, then results for a saxophone and a clarinet (Sec. 
m and finally the conclusions (Sec.|5]l. 



2 Background 



2.1 The TMM 



The transfer matrix method (TMM) provides an efficient means for calculating the input impedance of a hypothetical air 
column lisl ls. 121. With the TMM, a geometrical structure is approximated by a sequence of one-dimensional segments, such 
as cylinders, cones, and closed or open toneholes, and each segment is represented by a transfer matrix (TM) that relates its 
input to output frequency-domain quantities of pressure (P) and volume velocity ([/). The multiplication of these matrices 
yields a single matrix which must then be multiplied by an appropriate radiation impedance at its output. That is: 



p 



ZradU out 
Uout 



(1) 



where Zrad is the radiation impedance and Uout can be set arbitrarily. The input impedance is then calculated as Zj„ = 

Pin jUi-fi- 

The theoretical expression of the transfer matrix of a cylinder is: 



cosh {TL) Zo sinh {TL) 
smh{TL)/Zo cosh (TL) 



(2) 



where = pc/Tra?, F = jw/c + (1 + j)Q! is the complex propagation constant, p is the air density, c the speed of sound, and 
oj the angular frequency. Losses are represented by a, which depends on the radius a of the cylindrical pipe and varies with the 
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square root of the lossless wavenumber k ^ tu /c: 



a = 



(3) 



where CST is a constant that depends of the properties of air: 



CST = VW2(1 + (7 - l)/VPr), 



(4) 



£v = n/pcis the characteristic length of viscous effects, /i is the dynamic viscosity, Pr the Prandtl number and 7 the ratio of 
specific heats. 

The transfer matrix of a conical waveguide is (see Chaigne and Kergomard ||4|]): 



{0-2/ o-i) cos (kcL) — sin {kcL)/kxi jZc sin (kcL) 

Z^^ + (fc^2:ia;2)"^) sin (kjj) + [x^^ - s^^) cos (kJJ)ljk\ {ai/a2) cos [k^L) + sin(fc^L)/fcx2 



(5) 



where ai and 02 are the radii at the input and output planes, respectively, and xi and X2 are the distances between the apex of 
the cone and the input and output planes, Zc = pc/ {'Kaia2) and kc = —jT is the complex wavenumber In this case, losses are 
evaluated at the equivalent radius 141] : 



The transfer matrix of the tonehole is defined as: 



L 



ai 



1 



xi In (1 + L/xi) 



(6) 



- hole 



l + Za/2 Zail + Za/AZs) 
l/Z, 1 + Za/2 



(7) 



where Za is the series impedance and Zs the shunt impedance. These impedances have different values in the open and 



closed state. The calculation of these impedances has been the subject of many articles (Nederveen et al. llTll . Dubos et al 



list], Dalmont et al. 171], Lefebvre and Scavone ll4ll ) and the reader is referred to those papers for the appropriate formulas. 



2.2 The TMMI 

The radiation impedance of each tonehole on a woodwind instrument is influenced by the sound radiated from other holes. A 



method of solution to account for such interactions was proposed by Kergomard 11311 . It can be used for any bore shape by 
making use of the classical TMM with modifications for the matrices located between open toneholes. It gives identical results 
to the TMM if interactions are neglected (by specifying null mutual radiation impedances). That is, the geometry is discretized 
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Figure 1; Diagram of the flow contributions in Eq. (|9]l. 

identically, with both closed tonehole and open tonehole series impedance terms represented as in the TMM. 

We assume an instrument with openings (embouchure hole, toneholes, open end), where the indices of the openings 
range from n = 1 to A^. The pressure P„ inside the instrument at opening n is related to the acoustic flows C/„ radiating out of 
holes 71 by the following matrix relationship: 



P = ZV, (8) 

where we define the vector P of the pressures P„ and the vector U of the flow rates C/„ . Z is the radiation impedance matrix, 
which includes the effect of external interactions. The precise values of the different elements are difficult to determine. The 
self radiation impedances are the diagonal elements. The validity of this expression comes directly from the integral form 
of the Heknholtz equation if the Green function is chosen to satisfy the Neumann boundary conditions on the tube (see e.g. 
Eq. (7.1.17) in Morse and Ingard \v^, see also Leppington \15]). As a consequence, the equations used by Keefe |11| are 



erroneous (see Eqs. Ala to A2b). The content of this matrix is explained in Sec. 12.2.21 

An alternative equation relating the pressures and flows due to propagation inside the instrument can be derived by assuming 
conservation of acoustic flow at each opening n. As illustrated in Fig.[Tl the sum of the flow Un radiating out of the tonehole, 
the flow J7"^'** entering the tonehole section on the right, and the flow J/j'f ^* entering the tonehole section on the left is equal 
to the flow source [/* (which is discussed hereafter). This flow conservation equation can be written as: 

= Un + Ul^f' + C/^'s''*, (9) 



where we note that the flow on the left is defined in the reverse direction. 
This equation can be written in a matrix form as: 

Tie ft I -rrright 



U" = U + U"'--^' + (10) 



where U" is the flow source vector 
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The sum of the left and right internal flows for each section is related to the pressures: 

U'e/t _^U"S'" = YP, (11) 

where Y is the admittance matrix, which is described in Sec. 12.2.31 
By combining Eqs. (flOl) and (fTTT i , we obtain the solution: 

U = (I + YZ)-lU^ (12) 

where I is the identity matrix. 

2.2.1 Source Vector of Flow Rates U'' 

For such a calculation, the flow-source vector U'' needs to be known. Generally speaking, the reader can imagine a small 
loudspeaker located inside the pipe at the abscissa of each hole, providing a flow rate ?7*. Clearly the resonator of a musical 
instrument is passive and such sources do not exist. All transfer functions between two acoustic quantities at every point in the 
space of the passive system are fully determined and the unique problem is the choice of a reference. A solution is to use as a 
reference the flow rate on the left at the first tonehole of the instrument (from the part of the instrument that does not have any 
open holes), ~U^^^^ . In the absence of an active source, Eq. ( fTOb for this hole becomes: 

-U[^^' = Ui + ul'"'''. 

Thus, we can write Ul = —U^^^^ and replace Eq. ( fTol ) by the following: 

Ut^Ui + U{'^''\ (13) 

Otherwise = for n ^ 1. In this way, we can compute all quantities with respect to Uf, i.e. the ratio of all quantities to 
C/f . The input impedance can be easily deduced from the knowledge of Ul and Pi, the ratio U(/Pi being the input admittance 
y«p of the part of the system with open toneholes. Therefore the input impedance is classically computed by projecting the 
impedance l/y"P at the input of the instrument, or by using a transfer matrix relationship. Then all quantities can be calculated 
with respect to the input flow rate U^^^^* if necessary, where the index refers to quantities at the input plane of the system. 
For reed instruments, this quantity is related to the input pressure by a time-domain nonlinear characteristicQ 

'For flute-like instruments, it should be possible to choose the flow rate Ui exiting from the mouthpiece, which is the first open hole, as a source. However 
a complete nonlinear model needs to consider a pressure-difference source (i.e. a force source) near the edge, and it is necessary to add an equation in order to 
compute the flow rates radiating from the holes with respect to this source. 



5 



2.2.2 Radiation Impedance Matrix Z 

The self-radiation impedances Z„„ are approximately known (see e.g. Dalmont and Nederveen |16I]). On the other hand, the 



mutual radiation impedance Z^m (when n 7^ to) is, in its simplest form (see Pritchard ll9l Eq. (17)]) 



-jkd„ 



jkpc- 



-jkd„r, 



(14) 



where dnm is the distance between toneholes n and m. With this formula, it is assumed that a hole acts as a monopole. More 
closely spaced toneholes have a larger mutual radiation impedance. As the mutual impedance is complex, both reactive and 
dissipative effects are expected. 

Therefore, the impedance matrix Z is a full matrix. The mutual impedance may be neglected by using a diagonal matrix D 
with self impedance only, in which case the results are identical to those of the TMM. 

2.2.3 Admittance Matrix Y 

The propagation of planar sound waves between two toneholes 1 and 2 can be described by classical transfer matrices: 



Pn 
Jjright 



Cn D„ 



Pn+1 
left 



-u: 



n+l 



(15) 



where the transfer matrix is the multiplication of the transfer matrices of each segment located between the two open toneholes, 
including any closed toneholes. The series impedances of the open toneholes can be accounted for by including them in the 
transfer matrix, one-half on each side. 

This matrix can be written in the form of an admittance matrix: 



Tjright 
^ n 




rrleft 





Y Y 



Y., 



Y' 



Pn 
Pn+1 



(16) 



where the flow Ul^_ll is defined in the opposite direction. The parameters of this matrix are related to those of a classical 
transfer matrix: Yn = D„/i?„, Y^ = A„/i?„ and y^j.„ = — l/i3„, which assumes that A„_D„ — BnCn = 1, the condition for 
reciprocity. The right and left flows at one tonehole section n become: 



(17) 



and 



U, 



left 



^(1,71—1-^71—1 H~ ^fi—iPn 



(18) 
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Thus, Eq. (|9]l can be expanded to: 



^Un+ + + Yn)Pn + Y^^nPn+l. (19) 

The coefficients of this equation define the admittance matrix Y, which is tridiagonal. The first and last equations have to be 
modified because there is either no previous opening or no next opening. The last opening is located at the far end of the 
instrument, so that C/^'^'** = and Eq. ( fT9l l becomes simply: 

Uf, = UN + + YI^_^Pn, (20) 

where Un is the flow rate radiated at the end of the tube. 

For the first opening, we use Eq. ( fTSl l. where we can set JJf, the first entry of the flow source vector, to any value. Then, 
using Eq. ( fT2b . solving the problem gives the flow vector U. The pressure vector P can be calculated with Eq. 

2.3 Finite Element Calculations 

The evaluation of the input impedance of woodwind instruments using the FEM involves constructing a 3D model of the air 
column surrounded by a radiation sphere and the solution of the Helmholtz equation for a number of selected frequencies. The 
body of the instrument itself is considered to be rigid. The mesh occupies the volume inside and outside the instrument. Curved 
third-order Lagrange elements are used. 

The input impedance (or reflectance) is evaluated from the FEM solution by evaluating the relationship of pressure and 
volume flow (or traveling-wave components of pressure) at the input plane of the system. The surrounding spherical radiation 

n 

domain uses a second-order non-reflecting spherical boundary condition on its surface, as described by BayUss et al. 



Further discussion on this topic can be found in Tsynkov 



2l|] and Givoh and Neta 1911 



Thermoviscous boundary layer losses may be approximated with a special boundary condition such as presented by Cremer 
1 5] and, more recently, Bossart et al. |2] or Kampinga et al. ifioj . The boundary condition can be written as a specific acoustic 



admittance: 



1 



K-a« = -— = -^jkt, sin2 + (7 - 
p pc I 



(21) 



where w„ is the normal velocity on the boundary and 6 is the angle of incidence of the plane wave. The angle of incidence 
may be calculated from cos 6* = tt, • {)/| IwH, where the normal vector h is of unit length. This is solved iteratively. The lossless 
problem is solved first, then the admittance on the boundary is calculated from the normal velocity of the solution and the 
problem is solved again. This is repeated until convergence is found. 

The properties of air at 25 °C are used for all the simulation cases. See Causse et al. yfl for the equations used to calculate 
those values. 
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Figure 2: The mesh of a pipe with 10 open toneholes. 



The reflectance R{f) = p-/p+ (ratio of the reflected to incident pressure) is obtained from the simulation results. A 
cyUndrical segment is added before the input plane of the object under study. The pressures pa and pi, at two points on the 
centerline of this cylindrical segment, at distances a and b from the input plane, are extracted and the reflectance is calculated 
as: 

g-rb _ Hbae"^" 

^ = — Fa ° rfc ' ^^^^ 

where Hba = Pb/Pa is the transfer function between the two pressures and T is as previously defined. A singularity in 
this equation exists when the distance is half of the wavelength. The reduced impedance can then be calculated with Z = 

{1 + R)/{1-R). 

This method to calculate the reflectance was inspired by the two-microphones transfer function method of impedance 
measurement. It is worth mentioning that the impedance could also be calculated as Zin — Pin/pcvin, where the pressure and 
velocity are extracted directly at the input plane. When validating this approach using a cylindrical pipe, it was found that the 
results did not match theory as well as with the two-point method. 



3 Results for a Pipe with a Regular Array of Holes 

External tonehole interactions have been studied experimentally by Springer ll2oi during an internship in Le Mans, France. The 
experiment involved measuring the internal pressure at the position of the holes on a tube with an array of widely spaced tone- 
holes. The distances between the toneholes was much larger than what is found on woodwind instruments but the conclusion 
remains applicable to some extent. A cylindrical pipe of 4 meters length with an internal diameter of 15.3 mm was drilled with 
47 holes of 8.7 mm diameter regularly separated by 8 cm. The wall thickness was 3 mm. The internal pressure was measured 
at the positions of hole 1, 3 and 11. The transfer functions with respect to the pressure at the first tonehole were calculated. 
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Figure 3: Modulus and argument of the transfer function between the internal pressure at hole 3 and 1: experimental results 
(solid), theoretical calculation without external interactions (TMM, dotted) and with interactions (TMMI, dashed). 
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Figure 4: Modulus and argument of the transfer function between the internal pressure at hole 1 1 and 1 : experimental results 
(solid), theoretical calculation without external interaction (TMM, dotted) and with interaction (TMMI, dashed). 



The results are shown in Figs.[3]and|4]in comparison to theoretical calculations with and without interactions. 

For frequencies lower than the cutoff frequency of a tonehole lattice, the sound is exponentially attenuated inside the 
waveguide, whereas the external pressure is inversely proportional to distance. Therefore, the acoustic pressure coming from 
the outside of the toneholes located farther down an instrument becomes stronger than the pressure coming from inside the 
instrument. In Fig. [3] it appears that the effect of the external interactions is negligible for the 3rd tonehole because the 
pressure coming from inside remains important but, in Fig. |4] the internal pressure has sufficiently decayed at the 11th hole 
such that the external sound field dominates. The phase curve is very instructive: when interactions are ignored the phase 
shift is very small, indicating evanescent waves, while when interactions are taken into account, the phase variation is linear, 
indicating (spherical) traveling waves. 

For frequencies higher than the cutoff frequency of the tonehole lattice, the internal pressure is no longer exponentially 
attenuated and the effect of interactions is Umited to a smoothing of the response. 
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Figure 5: Magnitude of the reflectance (left) and equivalent length (right) of a pipe with 10 toneholes: experimental results 
(solid), FEM simulation results (squares), calculations with external interactions (TMMI, dotted) and calculations without 
interactions (TMM, dashed). 

In order to understand better the impact of the existence of external interactions on the playing characteristics of woodwind 
instruments, the influence of the external interactions on the input impedance (or reflectance) of a pipe with an array of closely 
spaced toneholes was studied. The pipe was 300mm in length, with a 12.7 mm diameter and 3.2 mm wall thickness, and was 
drilled with 10 holes of 9.5 mm diameter equally spaced by 15 mm. The reflectance of that pipe was obtained with the proposed 
calculation method and compared to simulation results with the FEM, to experimental measurements and to calculations with 
the classical TMM. The magnitude of the reflectance and the equivalent length are plotted on Fig.|5] 

An important observation is that the FEM results closely match the experimental measurement. This significantly increases 
our confidence in both the FEM and the experiment. For the equivalent length, the measurement error appears to be larger, 
particularly for the lower frequencies (below 1000 Hz). There also seems to be a slight systematic error of a few millimeters. 
The calculations with the proposed method to account for external interactions clearly gives better results than the TMM. The 
deep minima in the magnitude of the reflectance and the large increase in equivalent length in the higher frequencies completely 
disappear when interactions are included. The overall shape of the curves resemble the measured and simulated ones, even 
though some discrepancies remain. In the lower frequencies, the magnitude of the reflectance is reduced by the external 
interactions, which indicates a higher radiation efficiency. In the higher frequencies, the minima in the magnitude of the 
reflectance is not as low as in the measurement and is not located exactly at the same frequency. Small discrepancies also exist 
in the equivalent length, though they appear to be on the order of the measurement errors. In the lower frequencies, the external 
interactions increase the equivalent length slightly compared to predictions of the TMM. The toneholes on the pipe are located 
very close to each other, so that the evanescent modes excited near each discontinuity interact with those of adjacent toneholes, 
that is, the propagation of sound between toneholes is not planar, as assumed in the proposed method. This phenomena is one 
likely cause of the remaining discrepancies. Another is that the model of the mutual interaction assumes that each tonehole is 
a monopole. In spite of those simplifications, the proposed method gives improved results. Most of the discrepancies between 
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the classical TMM and the measurements are explained by the presence of external tonehole interactions. 

Generally speaking. Fig. |5] exhibits a major feature of stop bands: external interaction yields a significant reduction of 
oscillations with frequency, i.e. a reduction of the standing wave amplitude. This feature was stated by Kergomard Q. In 
Appendix A, a theoretical justification is given, allowing the following interpretation: 

• without interaction, there is reflection at the end, with standing waves inside the lattice; 

• without interaction, standing waves imply the existence of extrema of flow rate, the different holes radiating at different 
levels; 

• the holes radiating strongly have an important influence on the holes radiating weakly, thus there is a kind of equalization 
of radiation by the different holes, thus a diminution of the apparent standing wave ratio (SWR); 

• finally at the input of the lattice there is a diminution of the reflection coefficient. 

A consequence is the reduced height of the impedance peaks above the cutoff frequency and a reduction in the radiation 



directivity lobes in the backward direction 11 1311 . 



4 Results for a Saxophone and a Clarinet 

A precise computational FEM model of a complete music instrument is difficult to create and requires significant computation 
time to solve. Thus, the TMMI model can provide a faster and easier numerical technique which provides satisfactory results 
for real instruments with complicated geometry, despite the fact that the theoretical description of the toneholes with key pads 
is overly simplified. In general, we can at least expect that qualitative effects are well represented. Results of the TMM, TMMI, 
and measurements for an alto saxophone and a clarinet are presented in this section. 

4.1 Saxophone 

The input impedance of an alto saxophone was measured and compared with calculations using the TMMI and classical TMM 
methods. The instrument is a Sekner Super Action Series II, serial 438024. The imaginary part of the reflectance and the 
magnitude of the impedance for the first register and C#4 fingerings (respectively 247.0 and 329.6Hz) are shown in 
Figs.|6]and|7l 

For both fingerings, it appears that the TMMI better matches the experimental results than the TMM. The magnitude of the 
first three peaks of the impedance curve are reduced by the external interactions and closely match the experimental data. The 
behavior of the imaginary part of the reflectance is also better predicted by the TMMI. 
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Figure 6: Imaginary part of the reflectance (left) and magnitude of the impedance (right) of an alto saxophone with a Ab 
fingering: experimental results (solid) calculations with external interactions (dashed) and calculations with the TMM (dotted). 




Figure 7: Imaginary part of the reflectance (left) and magnitude of the impedance (right) of an alto saxophone with a CjJ 
fingering: experimental results (solid) calculations with external interactions (dashed) and calculations with the TMM (dotted). 

4.2 Clarinet 

Similarly to the case of the saxophone, the resonance frequencies of a Bb clarinet are predicted to be lower when external 
interactions are accounted for 

For fingerings where many toneholes are open, the lowering is on the order of 10-15 cents, slightly larger than for the 
saxophone. Of course, the lowest notes of the instrument, where only a few toneholes are open, are not much affected by this 
effect. 

For higher frequencies, the behavior of the instrument changes significantly. As an example, Fig.[8]displays the imaginary 
part of the reflectance and the magnitude of the impedance for the fingering A4 (391 Hz). The resonances correspond to the ze- 
ros of the imaginary part of the reflectance when the slope is negative. For the first two resonances, the external interaction only 
slightly shifts the frequencies to a lower value and the maximum of the impedance corresponds with the zeros of the imaginary 
part of the reflectance. There are three successive maxima between 1500 Hz and 2000 Hz. The frequencies and magnitudes 
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frequency (Hz) frequency (Hz) 

Figure 8: Imaginary part of the reflectance (left) and magnitude of the impedance (right) of a clarinet with a A fingering: 
experimental results (solid) calculations with external interactions (dashed) and calculations with the TMM (dotted). 

of these maxima are significantly modified by the external interaction. Notably, the first maxima no longer corresponds to a 
zero of the imaginary part of the reflectance so that only two resonances exists in this range compared with three when the 
interactions are accounted for. 



5 Conclusion 

The aim of this paper was the statement of the TMMI calculation method, comparisons between FEM, TMMl, and TMM 
methods, and preliminary comparisons with experimental data. Future work needs to be done for a systematic comparison 
between theory and experiment for the case of woodwind instruments. This is a long and delicate task because it requires 
very precise geometrical measurements, including bends and positions of the keys over the tone holes, as well as a precise 
theory. Concerning the TMMl, it would be useful to add internal interaction effects, at least for the determination of the 
resonance frequencies (notice that these effects do not yield radiation effects, thus dissipative effects, in comparison to the 
external interactions). 

The TMMl method provides a more accurate means for the calculation of the acoustics properties of woodwind instrument 
because it accounts for external interactions. The discrepancies between the TMMl and FEM are rather small and can probably 
be explained by several factors: first, the values of the radiation-matrix elements are roughly approximated; second, when 
several adjacent toneholes are closed, further improvement of the higher frequency modeling of a woodwind instrument would 
require internal coupling of higher-order modes to be accounted for 

In Appendix A, a modified version of the TMM is investigated. It is based upon a perturbation principle, and could be 
simpler than the TMMI. Unfortunately, the convergence is limited to pass bands (above the first cutoff). Therefore this method, 
used by Nederveen et al. 11171 p. 1 15], cannot be used for stop bands. 



Finally, we can summarize some effects of the interactions between holes as: 
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• The effect of tonehole interactions is generally more important when the toneholes are closer together. 

• The order of magnitude of the interaction effect seems to be of the same order for saxophones and clarinets. 

• At higher frequencies, the high-pass filtering behavior of the tonehole lattice allows more flow past the first tonehole and 
increases the effect of interactions, particularly near the cutoff frequency. 

• Above cutoff, the standing wave ratio is reduced by the effect of interactions. 
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Appendix A: Is it Possible to Compute the External Interaction by the Transfer 
Matrix Method? 

We consider the equation to be solved: 

U= (1+YZ^)"^U, (23) 

It is interesting to study if it is possible to solve this equation by perturbation, starting from the TMM method. A quite natural 
way to do this is to consider that the effect of the external interaction is weak, and to keep a calculation based upon transfer 
matrices. A first calculation is done without interaction, then the pressures are modified by calculating them with interactions 
taken into account. The perturbation calculation can be stopped here, but it is possible to iterate it: a new self-impedance is 
calculated as the ratio of the modified pressure to the unmodified flow rate, then the new flow rates can be calculated again 
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from the transfer matrix method with the modified values of the self-impedances. In practice, the iteration scheme is found 
to converge for almost all frequencies except low ones. This result is intuitive because in the stop band, the external sound 
pressure decreases proportionally to the inverse of the distance, while the internal pressure decreases exponentially, therefore 
the external interaction is more significant. 

It is possible to derive a criterion of convergence for the iteration procedure and, when it converges, it is possible to prove 
that the result is correct. This is done hereafter At each step n of the calculation, the transfer matrix method leads to the 
following relationship between the source Us, having a single non-zero element, Us{l), and the pressure and flow rate vectors, 
P(") andU^"); 

Us= U("^+YP("1 (24) 

The calculation is done by defining a diagonal matrix for the termination impedance of each hole (both the direct method and 
the transfer matrix method can be used): 

p(n)^ J3(n)-jj(n)_ (25) 

From the knowledge of the flow rate U^"^ the next value of the pressure p("+i) is deduced: 

P("+1^= Zi^U("). (26) 



The iteration equation is therefore found to be, with M = YZ^j: 

U("+i)=Us-MU("\ (27) 



The recurrence relationship leads to the following solution: 



^(-l)'M' 



.1=0 



(28) 



If the norm of the matrix M is less than unity, the recurrence converges to the solution ( 1231 ). the series corresponding to a 
Neumann series expansion. 

Now the starting point can be discussed. The first idea is to deduce the solution without interaction from the transfer matrix 
product : 

U(°)== (1 + YD)"^U„ (29) 



where D is the diagonal matrix of the self-impedances of Z/j. Another possibility is to start with U(''''= Us: this implies that 
the first pressure vector is built with the pressures created by a flow rate located at the first open hole. It can be concluded that 
the transfer matrix method can be used when the norm of the matrix YZj^ is less than unity. Because it can be verified that 
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this is not true in the stop band, the perturbation method unfortunately cannot be used in general. This confirms the intuition: 
looking at Fig. |5] it appears that the effect of external interaction can be very large in stop bands for holes very far apart from 
each other and the perturbation method cannot converge. 

Nevertheless, in pass bands we observe that convergence occurs rapidly when starting from Eq. ( |29] l. and even the first 
order, corresponding to a single perturbation step, is satisfactory. This observation thus justifies the reasoning given in Sec. [3] 
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